Morphology-Density Relation for Simulated Clusters of Galaxies 
in Cold Dark Matter Dominated Universes 
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ABSTRACT 

We present a model to investigate the formation and evolution of cluster galaxies using cos- 
mological high-resolution TV-body simulations. The TV-body simulations are used to construct 
merging history trees of dark halos. Gas cooling, star formation, supernova feedback and mergers 
of galaxies within dark halos are included by using simple prescriptions taken from semi-analytic 
models of galaxy formation. In this paper, we examine the merger-driven bulge formation model 
and represent the morphology-density relation of cluster galaxies at z = 0. We find that this 
morphological evolution model can explain the distribution of elliptical galaxies in the clusters 
well and cannot reproduce the distribution of SO galaxies. This result suggests that the elliptical 
galaxies are mainly formed by the major mergers, while, in the SO formation, the processes other 
than major mergers play important role. 

Subject headings: galaxies: clusters: general — galaxies: halos — galaxies: interactions 



1. INTRODUCTION 

It is well established that galaxy populations 
vary with the density of neighbouring galax- 
ies in clusters of galaxies (Dressier 1980). This 
morphology-density relation (hereafter MDR) in- 
dicates that dynamical processes that depend on 
environment of each galaxy mainly affect the fi- 
nal configuration of the stellar component. This 
impression seems to be supported by the Hubble 
Space Telescope images of clusters at intermediate 
redshifts, which show an abnormally high propor- 
tion of spiral and irregular types at z ~ 0.5 and 
increase of SO fraction toward the present time 
(Couch et al. 1994; Dressier et al. 1997). 

Some mechanisms that may transform one mor- 
phological type into another have been proposed, 
for example, ram-pressure stripping of interstel- 



lar medium of spirals by a intracluster medium 
(Gunn & Gott 1972; Fujita & Nagashima 1999), 
galaxy harassment by cumulative tidal interac- 
tions in clusters (Moore et al. 1996, 1998), and 
galaxy mergers (Toomre 1977). TV-body simula- 
tions have confirmed that merging disk galaxies 
produce galaxies resembling ellipticals as merger 
remnants (e.g. Barnes 1996), so that the galaxy 
merger is a favorite explanation for the pre- 
dominance of ellipticals and SOs in rich clus- 
ters (Kauffmann et al. 1993; Kauffmann 1996; 
Baugh et al. 1996). 

The purpose of our study is to check this merger 
hypotheses by investigating the MDR for the simu- 
lated clusters in the cosmological context. For this 
purpose, it is difficult in the present situation to 
use numerical simulations including both gravity 
and gas dynamics (e.g. Katz et al. 1992; Evrard et 
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al. 1994), because such simulations are expensive 
in CPU time and so only a limited range of pa- 
rameter space can be explored with an insufficient 
spatial resolution. 

Semi-analytic modeling of galaxy formation has 
already proved to be a powerful technique (e.g. 
Kauffmann et al. 1993; Cole et al. 1994). How- 
ever, we cannot identify the position of galaxies 
by such approaches, because these models follow 
the collapse and merging histories of dark ha- 
los by using a probabilistic method on the mass 
distribution based upon an extension of the Press- 
Schechter formalism (Press & Schechter 1972; 
Bond et al. 1991; Bower 1991; Laccy & Cole 1993) 

One approach to identify the position and ve- 
locity of each galaxy is to track the merging 
histories of dark halos by using TV-body sim- 
ulations and to combine them with the sim- 
ple prescriptions of the semi-analytic models 
(Roukcma et al. 1997; Kauffmann et al. 1999; 
Benson et al. 2000). Their schemes, however, do 
not deal with the substructures within dark halos. 
Since the dynamics within clusters may strongly 
affects the evolution of cluster galaxies (Okamoto 
& Habe 1999, 2000), we should use a new galaxy 
tracing method. 

In this paper, we adopt the galaxy tracing 
method provided by Okamoto & Habe (1999, 
2000), which enables us to trace the individ- 
ual galactic dark halos within dense environ- 
ments using high-resolution TV-body simulations, 
that is, we can obtain the three-dimensional dis- 
tribution of the galaxies within clusters. Our 
formula of modeling gas cooling, star forma- 
tion, supernova feedback, and galaxy merg- 
ers is directly taken from the previous semi- 
analytic work. We also adopt a merger-driven 
scheme for the production of galactic bulges 
and the way of the morphological classification 
based on bulge-to-disk ratios as earlier studies 
(Kauffmann et al. 1993; Baugh et al. 1996). 

In §2, we describe a brief outline of the galaxy 
formation model used here. We show the results 
about the morphology of the galaxies in §3. These 
results are discussed in §4. 

2. MODEL 

We exam the evolution of cluster galaxies in 
the standard cold dark matter (SCDM) universe 



(n = l,h = i/o/lOOkm s- 1 Mpc- 1 = 0.5, cr 8 = 
0.67) and the open CDM (OCDM) universe (fi = 
0.3, h = 0.7,(78 = 1-0). The baryon density is 
set to fl b = 0.1 and 0.03 for SCDM and OCDM, 
respectively. 

The outline of the procedures of galaxy for- 
mation is as follows. At first, the merging path 
of galactic halos are realized by the cosmological 
high- resolution TV-body simulations. Next, in each 
merging path, evolution of the baryonic compo- 
nent, namely, gas cooling, star formation, and su- 
pernova feedback, are calculated based on Kauff- 
mann et al. (1993) and Cole et al. (1994). We 
refer a system consisting of the stars and cooled 
galaxy. When two or more dark halos 
merge together, we estimate the merging time- 
scale based on the dynamical friction time-scale. 
When the merging of galaxies occurs, we change 
the morphology of the merger remnant by the type 
of the merger. Finally, we calculate the luminosity 
and color of each galaxy. Through the above pro- 
cedures, we obtain the morphological distribution 
of cluster galaxies. 

2.1. TV-Body Calculation and Merging of 
Dark Halos 

The merging histories of galactic dark halos are 
realized using the same method and simulation 
data by Okamoto & Habe (2000). Note that they 
trace halo-stripped galaxies as well as the galactic 
halos because halo disruption is probably due to 
lack of dissipative processes which are not included 
in the TV-body simulations (Summers et al. 1995). 
The merging histories are constructed with a 0.5 
Gyr time step in their paper, while here we adopt 
half of the time step for high redshifts (z ;> 2.0) 
at which the galactic halos form and merge vi- 
olently (Okamoto & Habe 2000). This improve- 
ment, however, hardly changes our results. 

2.2. Model of Galaxy Formation 

The following prescriptions are almost the same 
as the previous ordinary semi-analytic models. 

For simplicity, a dark halo is modeled as an 
isothermal sphere whose mass and radius are taken 
from TV-body data. The source of the diffuse gas 
in a halo is hot gas contained in its progenitor ha- 
los and in the accreting matter. The baryon frac- 
tion of the accreting matter is defined as /b = jf- . 
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We assume that the hot gas has the distribution 
that parallels to that of the dark matter with the 
virial temperature of the halo. When a galactic 
halo is tidally stripped, the hot gas in the halo is 
also stripped proportional to the amount of the 
stripped dark matter. 

The cooling time-scale, t coo i, is obtained as a 
function of the radius from the hot gas density 
profile, the temperature of the gas, and the cooling 
function A(T) as follows, 



7"cool ~ 



3 p(r) kT 



2 /im p n*A(T) ' 



(1) 



where /zm p is the mean molecular weight, n e (r) 
is the electron number density at the radius r, 
and k is the Boltzmann constant. Using the zero- 
metallicity and solar-metallicity cooling functions 
given by Sutherland & Dopita (1991), the cooling 
efficiency depending on the metallicity of the gas 
are calculated by interpolation and extrapolation. 

When a halo mass becomes more than dou- 
ble of the mass at the forming time, the dif- 
fuse hot gas contained in the halo is reheated 
by shock to the virial temperature of the halo 
(Somerville & Primack 1999). We refer this epoch 
as the reheating time of the halo. The cool- 
ing radius, r coo i, is defined as a radius at which 
the cooling time-scale, r coo i, equals to the elapsed 
time from the last reheating time of the halo, t. 
The hot gas that distributes between r coo i(t) and 
r C oo\(t + At) is cooled and added to the cold gas 
reservoir of the galaxy during the time-step, At. 

The star formation in disks is described by the 
following simple law, 



M* = 



M, 



cold 



r* 



T* = T. 



I T dyn 
r dyn 



(2) 



(3) 



where eq.(2) denotes the rate of stars newly 
formed, r° is the star formation time-scale of the 
galaxies which forms at z = 0, Td yn = is the 
dynamical time-scale of the galaxies, and in- 
dicates the dynamical time-scale of a halo which 
forms at z = 0. Here, we assume that the star 
formation time-scale is proportional to the dy- 
namical time-scale. We calculate r° yn according 
to the virial theorem and the spherical collapse 



model. The star formation time-scale at z = 0, 
t" , is a free parameter and we set to 5 Gyr in this 
paper. We have confirmed that this parameter 
docs not affect strongly on the morphology of the 
cluster galaxies. 

The feedback process by supernovae of massive 
stars has many uncertainties actually, therefore we 
adopt a simple description (Cole et al. 1994), 



AM r , 



^hcat 



v hoi 



M„At = At, (4) 



where V^ot and ahot are free parameters. When 
adopting a^ ot = 2, eq.(4) corresponds to the for- 
mula in Kauffmann et al. (1993), and Cole et al. 
(1994) used Qhot = 5.5 in their fiducial model. 

When two or more halos merge together, we 
identify a galaxy contained in the largest pro- 
genitor as the central galaxy of the new com- 
mon halo. Other galaxies are identified as satel- 
lites. These satellites merge with the central 
galaxy when the elapsed time from the last reheat- 
ing time exceeds the dynamical friction time-scale 
(Binney & Tremaine 1987), 



' mrg 



lnAGM sat 



(5) 



where rhaio and V c arc the radius and circular 
velocity of the new common halo, respectively, 
M sat is the total mass of a halo to which the 
satellite belonged as a central galaxy or the mass 
of stars and cold gas in the case that the satel- 
lite was a halo-stripped galaxy, and In A is the 
Coulomb logarithm, which is approximated as ~ 
ln(l + M,2 alo /M s 2 at ) (Somerville & Primack 1999). 
When the mass growth of the common halo satis- 
fies the reheating condition, the elapsed time from 
the last reheating time is recalculated because or- 
bits of satellites may be violently disturbed in such 
a case. 

When a satellite galaxy merges with a cen- 
tral galaxy and the mass ratio of the galaxy with 
smaller mass to that with larger mass is larger 
than /bulge, all stars of the satellite and disk stars 
of the central galaxy are incorporated with the 
bulge of the central galaxy. Then cold gases of 
both galaxies turns to stars in the bulge by star- 
burst. We adopt /bulge = 0.3 for SCDM accord- 
ing to the numerical simulations (Barnes 1996). 
Since early formation of the cluster prevent galaxy 
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mergers in OCDM, we use / bu igo = 0.2 for OCDM 
in order to reproduce the observed E+SO fraction 
roughly (see §3). 

Table 1: Feedback parameters for models 



Model Vhot (km/s) q hot 



A 


280 


2 


B 


140 


2 


C 


200 


5.5 







model B 




■ model C 




Sondoge et ol. (1985) 













-24 -22 -20 -18 -16 -24 -22 -20 -19 -16 



Chemical evolution is treated in almost the 
same way as described in Kauffmann & Chariot 
(1998). The instantaneous recycling approxima- 
tion is adopted. The amount of metals ejected 
from supernovae is characterized by y, which is 
heavy element yield for each generation of stars. 
The fraction / of the ejecta is ejected directly into 
the hot gas, and the rest is incorporated with the 
cold gas. We adopt y — 2Z© and / = 0.3, which 
are the same values as the strong feedback model 
in Kauffmann & Chariot (1999). The gas fraction 
returned by evolved stars, i?, is 0.25 in this paper. 
Simultaneously, the supernovae heat up the sur- 
rounding cold gas, then metals contained in the 
cold gas are also returned to the hot gas. The 
chemical evolution mainly affects on the colors of 
the galaxies. Therefore, the value of the yield is 
not important in this paper. The colors of the 
cluster galaxies will be discussed in our next pa- 
per. 

In order to compare our results with observa- 
tions directly, stellar population synthesis model 
must be considered. We use the model by Ko- 
dama & Arimoto (1997) with the Salpeter's IMF 
having a slope of 1.35 and mass range between 
O.IMq and 60M©. The range of stellar metallicity 
of simple stellar populations is 0.0001 ~ 0.05. 
The luminosity of disks and bulges is given by this 
model in each band according to the metalicities 
and ages of stars. 

2.3. Identification of morphology 

Morphology of each galaxy is determined by the 
S-band bulgc-to-disk luminosity idX\o(B / D) . In 
this paper, galaxies with B/D > 1.52 are iden- 
tified as ellipticals, 0.68 < B/D < 1.52 as SOs, 
and B/D < 0.68 as spirals, according to the re- 



Fig. 1. — Luminosity functions of cluster galaxies. 
The left and right panels show those in SCDM and 
OCDM, respectively, and the thick solid, dotted, 
and dashed lines indicate those of the model A, B, 
and C, respectively. The thin dot-dashed line is 
the luminosity function of the Virgo cluster galax- 
ies by Sandage et al. (1985). 

suits of Simicn & de Vaucouleurs (1986). It has 
been shown that this method for classification re- 
produces observations well by Kauffmann et al. 
(1993) and Baugh et al. (1996). 

2.4. Parameter Sets 

The feedback is key process which determines 
the feature of galaxies (e.g. Kauffmann & Chariot 
1998). Therefore, we use three type of feedback 
models in this paper (Table 1). The model A is a 
strong feedback model and B is a normal feedback 
model. In these two models ahot is set to 2. In 
the model C, we adopt ahot = 5.5 which is used in 
the fiducial model of Cole et al. (1994), while Vhot 
is equal to 200 km s~ x that is large compared to 
them. 

In Fig. 1, we show the -B-band luminosity 
functions of model galaxies. The thick solid, dot- 
ted, and dashed lines indicate the luminosity func- 
tions given by the models A, B, and C, respec- 
tively, and the thin dot-dashed line indicates the 
observed luminosity function of the Virgo cluster 
(Sandage et al. 1985). The model C reproduces 
the observed luminosity function well by flatten- 
ing the faint-end slope of the luminosity function. 
This is caused by strong feedback to galaxies with 
low circular velocities compared to the models A 
and B due to ahot = 5.5. 
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Fig. 2. — E + SO fractions of cluster galaxies as a 
function of absolute B magnitude for SCDM (left 
panel) and OCDM (right panel). The thick solid, 
dotted, and dashed lines indicate the model A, 
B, and C, respectively. The thin dash-dotted line 
shows the observational result for the Virgo cluster 
taken from Sandage et al. (1985) 

3. Morphology of the Cluster Galaxies 

3.1. Morphological Fraction 

In Fig. 2, the E+SO fractions of our clus- 
ter galaxies are shown as a function of absolute 
B-band magnitude. The observed curve for the 
Virgo cluster galaxies is taken from Sandage et 
al. (1985). All our models show the observational 
trend. It is, thus, found that the merger-induced 
bulge formation naturally increases the fraction of 
bulge-dominated galaxies toward a bright end. In 
the model C, the clusters do not have early type 
galaxies sufficiently. This is probably caused by 
the lack of SO galaxies as we show below. 

3.2. Morphology-Density Relation 

The MDRs for the simulated clusters are rep- 
resented in Fig. 3 using the same luminosity cut- 
off and definition of the local projected density as 
Dressier (1980), i.e. the local projected density is 
defined by nearest 10 neighbours having the lumi- 
nostiy My — 5 log h < —18.9. The projected den- 
sity is calculated in x—y, y—z, and z—x projections 
for each model. The morphological fractions of 
our models and Dressler's (1980) are represented 
by thick and thin lines, respectively. We attach 
the la error bars to the elliptical fractions accord- 
ing to the number of the galaxies in each density 
bin. 



In all models, the E fractions increase toward 
the high-density regions, and this is consistent 
with the observed trend. The SO fractions, how- 
ever, are much smaller than the observed frac- 
tion and do not represent the observed trend of 
the moderate increase toward the high-density re- 
gions. We then examine whether such deficiency 
of SO galaxies within the simulated clusters can be 
solved by adjusting the /bulge and the B/D clas- 
sification in the next subsection. 



SCDM OCDM 




Fig. 3. — Morphology-density relations for clus- 
ter galaxies. The left and right panels show the 
results for SCDM and OCDM, respectively. The 
thick lines are the morphological fraction of the 
thin lines are those of the observed cluster galax- 
ies by Dressier (1980). S, E, and SO fractions are 
represented by the solid, dotted, and dashed lines, 
respectively. We show the error bars for the ellip- 
tical fractions, which are la Poissonian uncertain- 
ties estimated from the number of galaxies in each 
bin. 
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Fig. 4. — Same as Fig. 2, but for the model Ci 
(thick solid) and C2 (thick dotted). For compari- 
son, the fraction in model C and the fraction of the 
Virgo clusters are represented by the thin dashed 
line and the thin dash-dotted lines, respectively. 

3.3. Dependence on Model Parameters 

Since there is considerable scatter in the rela- 
tionship between the bulge-to-disk ratio and Hub- 
ble T-type (Baugh ct al. 1996), we should examine 
the case in which we adopt different values of B /D 
for the classification of the galaxies. We also in- 
vestigate how the parameter, /bulge, affects on the 
MDRs. 

To study the influence of the choice of these pa- 
rameters, we use the model C in OCDM. By this 
model, the luminosity function and the distribu- 
tion of the elliptical galaxies are in good agree- 
ment with the observations, the E+SO fraction, 
however, is insufficient. 

First, we adapt the new B/D classification, 
that is, galaxies with B/D > 1.52 are identified 
as ellipticals, 0.25 < B/D < 1.52 as SOs, and 




o.oLuJ 1 1 1 1 1 l 

1 10 100 1000 

I Mpc" 2 



Fig. 5. — The MDR for the cluster in the model 
Ci. We show the la error bars for the SO fraction. 

B/D < 0.25 as spirals. We call this model the 
model Ci. In Fig. 4, wc plot the E+SO frac- 
tion of this model (solid line). The fraction is 
slightly increased, however, still below the obser- 
vation. The MDR of this model is presented in 
Fig. 5. Although the SO fraction is increased 
by adopting wider B/D range for SOs, the frac- 
tion is too small, especially in the high-density 
regions. Koopmann & Kenney (1998) indicated 
that galaxies with low central concentration (i.e. 
low B/D), with which they should be classified 
into spiral galaxies such as Sa galaxies in field en- 
vironment, were often identified as SO galaxies due 
to their low star formation rates. They suggested 
that misleading classification of the low central- 
concentration galaxies as SOs might account for a 
part of domination of SO galaxies in cluster envi- 
ronment. Our result shows that the distribution 
of the SOs is not reproduced by the merger-driven 
bulge formation even when we identify the galax- 
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ies with low B/D as SOs. 



4. DISCUSSION 



so 




Fig. 6. — The MDRs for the cluster in the model 
C2. In the left panel, we show the MDR adopt- 
ing the morphological classification mentioned in 
§2.3. In the right panel, we represnt the MDR in 
the case that we identify the galaxy with 0.68 < 
B/D < 4 as SOs. We show the la error bars for 
the SO fractions. 

Next, we set /bulge to 0.1 in order to increase 
the E+SO fraction. Note that this value of /bulge 
is extremely low. Indeed, a merger with the 
mass ratio of 0.1 should be considered as a minor 
merger (Walker et al. 1996). We call this model 
the model C2. The E+SO fraction of model C2 
with the B/D classification in §2.3 is shown in 
Fig. 4. The fraction is fairly increased than the 
model C and show better agreement with the ob- 
servation. In the left panel of the Fig. 6, we show 
the MDR of this model. It is found that both 
of E and SO fraction are increased, the shortage 
of SOs is, however, still significant. We then ex- 
amine the case in which we change the threshold 
value of the B/D between E and SO in order to 
increase the SO fraction. We identify the galax- 
ies with B/D > 4 as Es, and then the galaxies 
with 0.68 < B/D < 4 as SOs. The MDR in this 
case is indicated in the right panel of Fig. 6. The 
SO fraction becomes comparable with the observa- 
tion in the low-density regions, whether it hardly 
changes in the high-density regions. It means that 
the galaxies in very high-density regions mainly 
separate into two types, i.e. bulge dominated and 
disk dominated galaxies, under the merger hypoth- 
esis. 



In this paper, we represent a new method 
that combines high-resolution iV-body simula- 
tions with a semi-analytic galaxy formation model 
for cluster galaxies. The high-resolution simula- 
tions enable us to identify and trace galactic halos 
within clusters directly. Therefore, we can obtain 
the three-dimensional distribution of the cluster 
galaxies and also we can incorporate the dynami- 
cal processes, for example merging and tidal strip- 
ping of the galactic halos in the clusters, into our 
modeling of galaxy formation. For the first study 
of the cluster galaxies using this method, we have 
shown the MDRs for the simulated clusters of 
galaxies with two cosmological models based on 
the merger-driven bulge formation scenario. 

The model which reproduce the observed lumi- 
nosity function well (model C) is successful in ex- 
plaining the observed elliptical fraction as a func- 
tion of the local density. However, any models 
cannot reproduce the observed trend in the dis- 
tribution of the cluster SOs. In all models, the 
SO fractions are too small, especially in the high 
density regions. 

We also examine the case in which the SO galax- 
ies in the cluster environment have the different 
B/D range from the SO galaxies in the field, be- 
cause Koopmann & Kenney (1998) found that 
the galaxies with low central concentration, which 
should be identified as spiral galaxies in the field, 
were often classified into the SO type in clusters. 
We have confirmed that the SO fraction are in- 
creased only in the low density regions even if we 
adopt lower threshold value of B/D between S and 
SO. 

We then investigate how the parameter, /bulge, 
affects on the MDR. When we choose smaller 
/bulge, the early- type fraction is increased. In this 
case, however, the elliptical fraction is mainly in- 
creased, the shortage of SOs is then still signifi- 
cant. Even if we classify the low B/D ellipticals 
as SOs in the model with small /bulge, this change 
increases the SO fraction only in the low density 
regions again. 

Above results imply that, under the merger- 
driven bulge formation model, the galaxies sep- 
arate into two types, i.e. almost pure bulge and 
pure disk systems, in the high density regions in- 
dependent on the choice of the parameters and 
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cosmologies. 

The reason is probably considered as follows. 
In the model adopted here, a SO galaxy is only 
formed by the disk formation after the last ma- 
jor merger. The efficiency of merging is, however, 
high before cluster formation, and then it rapidly 
decreases due to the large internal velocity disper- 
sions of the clusters and the reduction of the size 
of tidally truncated halos as the mass of clusters 
grows. The accretion onto the galactic halos is also 
prevented by the strong tidal field after beginning 
of the cluster formation (Okamoto & Habe 1999). 
Hence, the cluster ellipticals formed by the merg- 
ers at high redshifts hardly change their morpholo- 
gies into SOs by minor mergers or gas accretion, 
which form additional stellar disks. 

Our results suggest that the merger hypothesis 
can give a good explanation of the distribution of 
cluster ellipticals. On the other hand, to repro- 
duce the observed SO fraction, the environmen- 
tal effects that change the morphology of spiral 
galaxies into SO galaxies directly should be con- 
sidered, for example the ram-pressure stripping 
(Gunn & Gott 1972; Fujita & Nagashima 1999), 
the galaxy harassment (Moore et al. 1996, 1998), 
and the minor mergers (Walker et al. 1996). 

The morphological classification based on B /D 
is the simplest way. There are various ways of 
classification other than the B/D classification. 
For instance, Caon & Einasto (1995) found an- 
other type of morphological segregation of early- 
type galaxies in the Virgo cluster, based on the 
configurations of their isophotes. Comparing local 
projected densities, galaxies with boxy isophotes 
are located in local density enhancements, while 
galaxies with disky isophotes lie in regions of lower 
local densities, independently of their classifica- 
tion as elliptical or lenticular. Although, it is 
interesting to investigate what physical processes 
are responsible for this segregation, more numeri- 
cal studies are needed to include such morpholog- 
ical classification in our model because we need 
to understand the dynamical states of the internal 
structures of galaxies to modeling it, and it is then 
left for further studies. 

In the next paper, we will investigate whole 
properties of cluster galaxies in detail, i.e. col- 
ors, metalicities, velocities, and so on, using the 
method presented in this paper. We will also dis- 
cuss some environmental effects. 
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